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We consider the Emery model of a Cu-O plane of the high temperature superconductors. We 
show that in a strong-coupling limit, with strong Coulomb repulsions between electrons on nearest- 
neighbor O sites, the electron-dynamics is strictly one dimensional, and consequently a number of 
asymptotically exact results can be obtained concerning the electronic structure. In particular, we 
show that a nematic phase, which spontaneously breaks the point- group symmetry of the square 
lattice, is stable at low enough temperatures and strong enough coupling. 

PACS numbers: 



Immediately following the discovery of high temper- 
ature superconductivity in the cuprates, it was realized 
P, that the novel physics of these materials is domi- 
nated by the strong, short-range repulsion between elec- 
trons. However, there has been considerable debate over 
what is the simplest "paradigmatic" model that captures 
the essential physics of the problem. Despite the fact, 
pointed out early on by Emery .2| and by Varma and 
coworkers (J] , that the minimal model which captures the 
essential local chemistry of the doped copper-oxide planes 
is the three-band copper-oxide or Emery model (defined 
below) , it has generally been the accepted practice among 
theoreticians to, instead, consider the single-band Hub- 
bard or t-J model - certainly reasonable models for study- 
ing the interplay between the localized quantum antifer- 
romagnetism of the undoped system and the charge de- 
localization produced by doping. Moreover, since none 
of these strongly interacting models can be solved in two 
dimensions (2D), any theoretical results that can be es- 
tablished with an acceptable degree of rigor can shed light 
on the observed physics of the actual materials. 

In this paper we show that there exists a limit (which is 
not wildly unphysical) of the Emery model about which 
a number of exact statements are possible. Specifically, 
despite the fact that the model itself possesses the sym- 
metries of the square lattice, the electron dynamics is 
quasi-one-dimensional in this limit. It is also possible to 
establish the existence of various electronic liquid crys- 
talhne phases 0, including especially an Ising- nematic 
phase I55 which spontaneously breaks the four-fold rota- 
tional symmetry of the underlying lattice. 

I. THE MODEL 

We consider a model defined on the copper-oxide lat- 
tice, shown in Fig. ^ the corresponding Hamiltonian 
operator is written explicitly in Appendix A. The copper 
sites define a simple square lattice at lattice positions R 
with lattice constant, a, while the oxygen sites sit at the 
center of the nearest-neighbor bonds on this lattice, and 
so define a second square lattice, rotated by 45° relative 



to the Cu lattice, with lattice positions R + {a/2)ex and 
R + {a/2)ey, which we will call site {R,x) and {R,y), re- 
spectively. The vacuum is defined as the state in which 
all the O p orbitals and Cu d-orbitals are full. The rele- 
vant Fock space is constructed by adding holes (removing 
electrons) from the Cu 3dx2-y2 and O 2po- orbitals {i.e. 
the 2i>,j; orbital associated with the oxygens at R, x and 
the 2pj^ for the oxygens at R, y.) The corresponding hole 
creation operators are dt, and pt, 
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FIG. 1: Schematic representation of the Cu — O lattice. The 
full circles represent Cu sites and the open circles are O sites. 
The various terms in the Emery-model Hamiltonian are rep- 
resented, as discussed in the text. 

The various interactions in the model (also shown in 
Fig. nj are defined as follows: The repulsion between 
two holes on the same site is Ud and Up, respectively, 
for a copper and oxygen site, while the repulsion be- 
tween two holes on an adjacent copper and oxygen or 
a nearest- neighbor pair of oxygens are Vpd and Vpp. All 
further neighbor interactions are neglected. The hop- 
ping matrix elements which transfer a hole between a 
nearest-neighbor O - Cu pair is tpd = t, while that be- 
tween nearest- neighbor O's is tpp, and the difference be- 
tween the energy of an electron on an O and Cu site is 
e > 0. The signs of the various hopping matrix elements 
are determined by the symmetry of the relevant d and 
p orbitals. However, a simple gauge transformation with 
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wave- vector tt = (vr/a) (1, 1), changes the signs so that 
all the relevant hopping matrix elements are positive. 

The insulating parent state of the undoped cuprates 
has one hole per unit cell which, because e > 0, live pref- 
erentially on the Cu sites. Additional doped holes, whose 
concentration per unit cell we denote x, go preferentially 
on O sites because Ud ^ e. 
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FIG. 2: Schematic representation of various states discussed 
in the text. The dark circles are Cu and O sites occupied by 
a hole; the open circles are sites not occupied by holes. 



II. THE STRONG COUPLING LIMIT 

We start by defining the "strong-coupling limit" of this 
model, in which the interaction strengths are large com- 
pared to the one-electron energies. Here we discuss the 
salient features of this regime. In Appendix IXI we present 
details and prove that the statements we make here are 
asymptotically exact in this limit. 

By the strong coupling limit we formally mean that 
we consider the model in the limit U/t oo, where all 
the interactions, Up,Ud-,Vpd,Vpp ^ U and the hopping 
matrix elements, tpd = t and tpp ^ t. Even in this limit, 
the physics depends on the finite ratios of the various in- 
teraction strengths. In particular, we will always assume 
that the following inequalities are satisfied: f/d > e > 0, 
and Ud > Up > Vpd > Vpp > 0, consistent with chemical 
intuition. (Somewhat more restrictive inequalities must 
be assumed in order to prove all the stated results, as 
is discussed explicitly in Appendix A.) Finally, since the 
hopping matrix elements depend exponentially on sepa- 



ration, we set tpp/t — > 0. This final assumption may not 
be well satisfied in the actual materials, where cluster 
calculations suggest that tpp/t ^ 1/3. We will study 
this model as a function of x and for arbitrary ratio of 
e/t. 



For the undoped system, x = 0, the ground-state has 
zero energy and is 2^ fold degenerate, with one hole on 
each copper. Of course, this degeneracy is resolved for fi- 
nite interaction strengths when antiferromagnetic super- 
exchange interactions, with J « 8t^ /UpV^^, are included. 
However, in the strong-coupling limit, this (and most of 
the other spin-physics we will encounter) involves energy 
scales that vanish as t/U and t/V 0; we will therefore 
ignore this physics at first, and then return to it when we 
consider "t — J-like" physics that arises from low order 
corrections to the strong-coupling limit. 

Neglecting the spin degeneracy, the first excited state 
is an exciton, shown in Fig. with a large energy 
Eex = Vpd + e + 0{t), and so can be ignored at low 
energies and temperatures. 

Now, consider one additional doped hole. Since U » 
V, adding one hole means increasing the number of oc- 
cupied sites by 1, and so necessarily costs a minimum 
energy of fi = 2Vpd + £■ Some possible representative 
states are shown in Figs. Et-d: 

• I^D shows the bare hole-state, which to zeroth order 
in t has energy fj,, 

• |2h shows a hole-exciton bound state, with zeroth 
order energy fi + e. 

• |21i shows a hole broken into two charge e/2 solitons 
separated by L = 4 sites, with zeroth order energy 

-I- Le. 

• shows a bent hole-exciton bound-state which is 
the lowest energy state which involves a distur- 
bance outside of this row and which has zeroth or- 
der energy -t- Vpp + e, and so can be neglected in 
the strong coupling limit j3|. 

The most salient point to notice is that if the doped 
hole is added to an oxygen on a given row, all states with 
energy near /i involve disturbances which are confined to 
the same row. Any state which involves a non-colinear 
disturbance, such as the bent-hole exciton in costs 
infinite energy in the strong coupling limit. 

In this limit, therefore, the number of holes on each 
row and each column of the lattice are separately con- 
served quantities, and the charge dynamics is purely one- 
dimcnsional. More precisely, in Appendix ^ we show 
that, to leading order in the strong coupling expansion 
and for tpp — 0, each row p has a conserved quantity Xp, 
and each column q a conserved quantity Yg, which quali- 
tatively correspond the number of hole quasiparticles on 
that row or column. Indeed, doped holes on distinct par- 
allel rows do not interact with each other. However holes 
on rows and columns do interact with each other (where 
they meet). We will show below that these interactions 
play a crucial role. Of course, when we back off from the 
strong coupling limit, or if we include a small but non- 
zero tpp, small effective interactions which violate these 
conditions will be generated. 
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III. THE ID DYNAMICS 

Consider a system in which we add a fixed number 
of doped holes to one and only one row of the lattice. 
Because the number of electrons in each row is con- 
served, none of these holes can leak out onto other rows 
or columns, which thus remain undoped. Indeed, the 
electron dynamics along this row is exactly equivalent to 
those of the ID Cu-0 model, which was analyzed previ- 
ously in Ref. 0: 

Hrow = ^[c] ,^Cj+i,o- + h.c] + ^ ejflj 

^[t/j-rij- tfij,x + Vpdfijfij+i] (1) 

3 

where even number sites are Cu sites and odd numbered 
sites are O, e2j = 0, e2j+i = e, U2] = Ud-, U2J+1 = Up, 

fij^a- ~ Cj^Cj^a, and hj = J2cr^3,<^- ^^^^ strong cou- 
pling Uj 00 limit, the charge degrees of freedom can 
be treated as spinless fermions 8], with effective Hamil- 
tonian 

He = -t^[c]cj+i + h.c] + '^[ejnj + Vpdnjfi-j+i], (2) 
i j 

while again the dynamics of the spin-degrees of freedom 
are obtained only when corrections to the strong coupling 
limit of order xt^ jUj are included. The density of spinless 
fermions per site is simply (2^^^)^^ fij = 1 -\- x. 

Manifestly, for a; = 0, the system is insulating, with 
one hole on each copper site and a charge gap Ac = 
2^(i + 0{t). For small, positive x, the doped holes are 
dilute and can be treated within the context of an effec- 
tive mass approximation, as free spinless fermions with 
creation energy Ac and effective mass, m* . In particular, 
the ground state energy per site (with W = h^ii^ /Qm*)\s 

E ^ Eq^ I^cx + Wx^ + 0{x^). (3) 

Both Ac and W are continuous functions oi e/t: Ac = 
2Vpd + tFi^{e/t) and W = {8tt^ / 3)10^ Fw (e/t). For e > t, 
Fa = 1 + 0(i/e)^ and W = 8(t/e)[l + 0{t/ef]. In the 
opposite limit, t/t — > 0^, the fermions fractionalize to 
form twice as many charge e/2 solitonic Fermions. How- 
ever, the ground state energy has the same x dependence, 
but with Fa = -4(i/e)[l 0(e/i)] and W = l + 0{e/t). 
An important qualitative point to recognize here is that 
the Fermi pressure is a decreasing function oi e/t which 
vanishes as e/t 00. Various correlation functions can 
be accurately estimated, as well, from the well known 
theory of the ID Luttinger liquid. 

For X = 1, the system is again insulating, with one 
hole on each site. Expanding about this limit, for 1 ~ x 
small, yields a result similar to those obtained for small x. 
Other interesting states occur in the vicinity of various 
commensurate values of x. For instance, a; = 1/2 cor- 
responds to commensurability 2 in the spinless fermion 



problem, where an incompressible CDW state is the 
ground state of this strongly interacting problem. 

Indeed, given the large number of exact, or well con- 
trolled approximate results that can be obtained for the 
IDEG, and the ease with which quantum Monte-Carlo 
simulations can be employed to flesh out the analytic re- 
sults quantitatively, we consider the problem of a single, 
Cu-0 row to be a solved problem. This also means that 
a large number of "fully nematic" states can also be com- 
pletely characterized. A fully nematic state is defined to 
be one in which which doped holes are placed only on 
rows (or only on columns). Since holes on neighboring 
rows do not interact at all (unless we were to add longer 
range interactions to the model), the dynamics of the 
holes on each row are determined by the same ID Hamil- 
tonian we have just analyzed. This does not constitute 
a complete solution of the problem, since states in which 
some doped holes lie on rows, and others on columns are 
still complicated, and require additional analysis to char- 
acterize. However, we will show below that, under many 
circumstances, the ground state is fully nematic. 



IV. THE NEMATIC PHASE 

We now move from the analysis of the hole dynamics 
along a single row or column, to study the phases of the 
full two-dimensional model. 



A. The nematic insulator, x = 1 

At X = 1 there are two holes per unit cell, and in 
the ground state each Cu and O site along each row is 
occupied by exactly one hole, while the O sites along 
columns are empty. The energy per unit cell of this state 
is 2Vpd + e. There is a second, degenerate, ground state 
obtained by exchanging rows and columns. Each ground 
state spontaneously breaks the the 90° rotational invari- 
ance of the square lattice but is translationally invari- 
ant since all unit cells are equivalent. Furthermore these 
strong coupling ground states at x = 1 have a charge 
gap. Hence this phase is a nematic insulator. 



B. The "classical limit," e/t 00 

In this limit (still with Ud > e) the charge degrees 
of freedom define a classical lattice gas which can be 
mapped precisely to an antiferromagnetic Ising model on 
the oxygen lattice with exchange coupling Vpp/A. Here 
spin- up indicates an occupied state and spin down an un- 
occupied state. Under this mapping, the magnetization 
of the Ising model is m = 1 — x and the Neel state is the 
insulating nematic state. 

The phase diagram of this model, shown in Fig. |3| is 
well known 9]. For a; < 1, there is a continuous finite 
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temperature transition from a high temperature disor- 
dered phase to the low temperature Neel phase. How- 
ever, at a critical x = Xc, there is a tricritical point, such 
that for X < Xc the transition is discontinuous. At low 
temperatures, for any < x < 1, there is two-phase coex- 
istence between a ferromagnetic (a; = Mott-insulating) 
and a Neel (a; = 1 insulating nematic) phase. 

At elevated temperatures, this classical phase diagram 
is relatively insensitive to the addition of perturbations 
to the Hamiltonian. Specifically, even for t/t > 0, the 
phase diagram is hardly altered by quantum effects so 
long as r ^ t/^Jt'^ + e^. However, at low temperatures, 
additional longer range interactions, either added explic- 
itly to the model or induced by quantum fluctuations, 
can affect the nature of the stable phases substantially. 
For instance, even at the classical level, including the ef- 
fect of weak Coulomb repulsion between holes on second- 
neighbor O orbitals will stabilize an electronic crystalline 
phase in a narrow range of x near 1/2 and at low enough 
temperatures. Here the doped holes form a period 2 
(Wigner crystal) density wave along each row, while the 
O's along vertical bonds remain undoped. More gener- 
ally, at low T and intermediate x, the phase diagram is 
complex and dependent on details. However, for x near 
or 1, we will see that quantum effects generically stabilize 
homogeneous quantum nematic phases, shown as shaded 
areas in Fig. 13 




a; 1 



FIG. 3: Schematic phase diagram as a function of temper- 
ature and doped-hole concentration in the strong coupling 
limit (with magnetic ordering suppressed) . Only the classical 
phase boundaries are shown; thick (thin) lines: discontinu- 
ous (continuous) transitions; circle: tricritical point; shaded 
regions: nematic phases stabilized at t/e > 0. 



C. Quantum effects, t/e > 0, for x and 1 — a; <C 1 

At low temperatures, T <C , quantum effects 

are important, even for small t/e. We start by addressing 



the nature of the low temperature phase for a; ^ 1. 

Since the doped-holes are fermions, one would gener- 
ally expect an associated Fermi pressure which tends to 
favor a uniform state. As shown in Eq.©, the ID dy- 
namics of the doped holes imply that the Fermi pressure 
^ Wx^, rather than the stronger Wx^ dependence of a 
2D Fermi gas. However, absent any attractive effective 
interactions between doped holes the Fermi pressure is 
always sufficient to stabilize a uniform phase at small x. 
(Effective attractions can arise from magnetic fluctua- 
tions but are negligible at strong coupling.) 

More specifically, an upper-bound to the ground state 
energy per site, E{x), at all x is given by the energy of 
the fully nematic state, Enem{x), i.e. the ground state 
energy of the ID Emery model computed from Eq.QJ 
or Eq.(|2Jl. It can easily be imagined that there is an 
isotropic state (or a less nematic nematic state) that has 
lower energy. However, while the properties of such an 
isotropic state cannot be computed exactly, a simple es- 
timate shows that for x ^ 1, the ground state is always 
nematic. 

To see this, we compare the energy of the nematic 
phase with that of an isotropic version of this phase. If 
we ignore the interactions between doped holes on rows 
and columns, the isotropic phase would have lower en- 
ergy, since there are twice as many ID systems each 
with 1/2 the density of doped holes, resulting in a fac- 
tor of 4 reduction in the Fermi pressure. Now, for x 
small enough, the contributions to the ground state en- 
ergy from the Vpp coupling between holes on crossing 
Cu-0 rows and columns is a is regular function of x, 
free of infrared divergences. In particular, when two 
doped holes approach each other at the intersection of 
a row and column, the repel each other strongly, and 
the probability of such an interaction is proportional to 
x^. Thus, the energy per site of the nematic state is 
-E-ncm = Eq + Ac x + Wx^ + . . ., whilc that of the isotropic 
state is Ei^o = -Bo + Ac x+{l/4)Wx^ + V''^x^ + ... where 
V'^^ is an effective repulsion between holes on intersect- 
ing rows and columns. Manifestly, at small enough x, the 
nematic state has lower energy. 

The Fermi pressure similarly stabilizes the uniform 
phase for x < 1. Here the nematic character of the re- 
sulting uniform phase is considerably more obvious: the 
nematic phase consists of an array of Luttinger liquids 
(one per row) with an effective Luttinger charge param- 
eter Kc and an effective charge velocity Vq. 



V. CORRECTIONS TO THE STRONG 
COUPLING LIMIT 

Low order corrections to the strong coupling limit re- 
solve the spin-degeneracies we have neglected, and pro- 
duce other important changes in the physics, which we 
will discuss elsewhere; here we comment on a few salient 
features. 
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A. Magnetic Interactions 

In the undoped Mott insulator at a; = 0, the most obvi- 
ous induced interaction is the super-exchange interac- 
tion between nearest-neighbor Cu spins, which to leading 
(4*^) order in t is J = U^/[{Vpd + efiUp + 2e)]{l + {Up + 
2€)/{2Ud)}- However, unlike the one-band model, under 
some circumstances, for the strong coupling limit of the 
Emery model other higher order interactions can be com- 
parable in magnitude to this interaction. Thus, there is 
a 4-spin ring exchange interaction on a plaquette gener- 
ated at 8*'' order, which does not vanish for J7 — > oo. 
(For Ud = Ud = oo and Vpd = Vpp — V e, we get 
J4 = {23/2)t^ /V"^ while the leading order contribution 
to J is J ^ t*tpp/V*.) Thus, whereas the Mott insu- 
lating ground state at x = of the single band model is 
inevitably magnetically ordered, the Emery model has, in 
addition, quantum disordered (likely dimerized) phases. 

B. 2D Charge Dynamics 

For any of the conducting phases discussed above, the 
corrections to strong coupling not only resolve the spin 
degeneracies, but also lead to important changes in the 
charge dynamics at asymptotically low energies. In par- 
ticular, under most circumstances, we expect that the 
peculiar non-Fermi liquid behavior resulting from the 
strictly ID dynamics of the strong coupling limit will be 
destroyed at vanishing temperatures by induced interac- 
tions (either proportional to tpp or /Vpp) which permit 
holes to hop from a row to a column. These couplings 
will either lead to a crossover to Fermi liquid behavior 
at low temperatures, or to a broken symmetry ground 
state. However, so long as these interactions are weak, 
they cannot restore the point group symmetry, so the ne- 
matic character of the resulting states should be robust 
[Tol Moreover, as is characteristic of quasi ID sys- 
tems, the non-Fcrmi liquid character will still be manifest 
at non-zero temperatures, energies, or wave- vectors. 

C. Coulomb interactions 



electrons {x < 0) remove holes from the Cu sites; these 
"doped electrons" have no local tendency to be dynami- 
cally confined to rows or columns. 



VI. CONCLUSIONS 

It has long been accepted that the Emery model pro- 
vides a reasonable description of the relevant electronic 
degrees of freedom in the cuprates. In this paper, we 
have obtained theoretically well controlled results in a 
strong coupling regime of this model. Although cluster 
calculations suggest that this limit may not be entirely 
appropriate in the real materials, the insights obtained 
here may nevertheless capture important features of the 
physics. 

Perhaps the most salient feature of the results obtained 
here is the existence of a strongly nematic phase in a sig- 
nificant portion of the phase diagram. This contrasts 
with the behavior of the same model in the weak cou- 
pling limit, where it behaves in similar fashion to the 
single band Hubbard model in which such a phase, if it 
occurs at all, is confined to special fillings associated with 
the proximity to van Hove singularities [TH Il2i | . Exper- 
imental evidence of the existence of a nematic phase 
in the cuprates was recently reviewed in Ref. |l3| . The 
other feature of the phase diagram is the existence of 
a large region of two-phase coexistence - phase separa- 
tion. This may be the simplest example of a generic 
tendency of highly correlated systems to form inhomoge- 
neous states [lj|. 
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The fact that the model considered has regions of two- 
phase coexistence, and others where various susceptibili- 
ties (including the compressibility) are large, means that 
the low temperature physics can be strongly modified 
by the effect of even weak additional interactions. Of 
these, the most obvious is the long-range Coulomb in- 
teraction, which always frustrates phase separation and 
instead results in various locally inhomogeneous phases 
such as stripe and bubble phases. 

D. Particle-Hole Asymmetry 

At least at strong coupling, the Mott insulating state 
at a; = is strongly particle- hole asymmetric: added 



APPENDIX A: DETAILS OF THE DERIVATIONS 

In this Appendix, we formalize the statements concern- 
ing the strong coupling limit. The Cu sites are labeled 
by the Bravais lattice vectors, R, and the hole density 
on the Cu is n{R). The hole density on the two O sites 
in the same unit cell are labeled n±x{R) and n±y{R), 
respectively. (Notice that n-.x{R) = n^^R — ax), and 
n^y{R) = ny{R — ay).) In the strong coupling limit, the 
Hamiltonian for the model described in Fig. ^ can be 
written as the sum of two terms, H = Hq + Hi, with 
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and 



(Al) 



R.cr 



where tr = ± is the spin label. As discussed in the text, 
we set direct 0-0 hopping amplitude tpp to zero. 

To begin with we consider just the effect of the unper- 
turbed Hamiltonian Hq. To make the proofs simpler, we 
will consider the case in which the following inequalities 
are satisfied: 

Ud>2Vpd, Up>2Vpd, Vpp>^. (A3) 

In this case, Hq is positive semidcfinitc, as we will now 
demonstrate. It is apparent from the structure of the 
Cu-0 lattice, that for every Cu site R it is possible to 
define four triangles, located respectively NE, NW, SW 
and SE of the Cu site, and each having the Cu site and 



two adjacent O sites for its vertices (see Fig. We 
will label a triangle {R, s, s'} according to its Cu vertex 
R and by a pair of labels s = ± and s' = ±, where 
{s, s'} = {+, +} corresponds to the triangle NE of R, 
{s,s'} = {-I-,—} to the triangle NW, etc. We denote by 
nss'{R) the total hole occupancy of triangle {R, s, s'}: 

Uss'iR) ^ niR) + ns^R) + ns'y{R). (A4) 

By using this notation we can equivalently write Hq in 
terms of the hole occupancy of each triangle and of the 
occupancy of the Cu and O sites: 



Ho = {^-Vr_ 



pd 



J2n{R) {n{R) 



- 1 



P T/ 



E {MR) (MR) - l) + MR) [ny{R) - l) } 



^ ^ (n,AR)-l) (nsAR)-2 



R,s,. 



Vpd 



Vpp~-YjYl {nx{R)ny{R) + ny{R)n^^{R) + n^^{R)n^y{R) + n^y{R)n^{R)^ 

R 



(A5) 



Since each operator appearing in this expression is indi- 
vidually positive semidefinite, so is Ho, so long as all the 
inequalities of Eq. ljA3|l are satisfied. 

(Strictly, in the strong coupling limit, violating the in- 
equalities of Eq. IA3I can lead to major restructuring of 
the ground state. For instance, for Vpp < Vpd/2, the sys- 
tem phase will separate for any x > 0, albeit if t is not in- 
finitesimal, quantum effects may lead to inhomogeneous 
ordered states. Indeed, close to the "fully frustrated" 
point Vpp = Vpd/ 2, interesting forms of quantum order- 



from-disorder effects can arise. We will not consider these 
interesting issues further in this paper. ) 

Thus, provided the inequalities of Eq. (jA3p are satis- 
fied, the Hilbert space of (generally degenerate) zero en- 
ergy states consists of the set of configurations in which, 
a) the Cu and the O sites are either empty or singly occu- 
pied (by holes) , b) each triangle is either singly or doubly 
occupied (also by holes), and c) nearest-neighboring O 
sites are not simultaneously occupied by holes. All other 
states are separated from these zero energy states by a 
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finite energy gap. This is the low energy Hilbert space of 
states that we will consider here. 

Hence, so long as 1 > a; > 0, there exist zero energy 
ground states and an extensive ground-state degeneracy. 
Here, the density of "doped holes," x, is defined in terms 
of the total hole density per site, 



02 



— Q O P 



a) 



which is, of course, a conserved quantity. {N is the total 
number of Cu sites on this lattice.) 

To analyze the strong coupling limit to lowest order 
in perturbation theory, we confine our attention to the 
zero energy subspace of the full Hilbert space. Thus, 
the effective Hamiltonian, Hcs, is obtained by simply 
taking matrix elements of the perturbing Hamiltonian, 
Hi, defined in Eq. (|A2|I . between states in this subspace. 
Higher order terms in the perturbative expansion can be 
obtained, as is done in deriving the t — J model from 
the large U limit of the Hubbard model, by including the 
effects of virtual transitions to the finite energy states of 
the unperturbed Hamiltonian. For the present purposes, 
we discuss only the first order problem, i.e. 



i 



H, 



off 



PoHiPo 



(A7) 



where Pq is the projection operator onto the subspace of 
zero energy eigenstates of Hq. (See Fig. 

For a; = 0, there is a three-fold orbital degeneracy of 
the unperturbed ground state, in addition to the 2^- 
fold spin degeneracy (which is only lifted, as discussed 
in the text, when high order superexchange processes are 
included). Any e > eliminates this orbital degener- 
acy, and uniquely choses the ground-state with one hole 
on each Cu site; there are no non-zero matrix elements 
of the hopping term between states in the ground-state 
manifold. 

For X — 1, there is a two- fold orbital degeneracy of the 
ground state, as already discussed in the text. As the 
two states involved are related by a symmetry operation 
of the Hamiltonian (rotation by 7r/4) this degeneracy is 
not lifted in any order of perturbation theory - this is the 
nematic insulating phase discussed in the text. 

For e ^ t, we perform a second perturbation expan- 
sion, which we refer to in the text as the "classical" limit: 
To first order in e, the ground state degeneracy is reduced 
to the subset of states which have precisely xN holes on 
O sites, i.e. the fewest number possible. Since the term 
proportional to t in Hi changes the number of holes on O 
sites, there are no matrix elements between these states 
to fist order in t. However, to order t^/e, it is possible 
for holes to move without violating this constraint. 

For intermediate values of x and t/e, the effective 
Hamiltonian is generally fairly complicated and we have 
not obtained a general solution. However, we will now 
prove that under the dynamics of H^jf, there are con- 
served quantities, Xp and Yp, corresponding to the num- 
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FIG. 4: Conservation of Xp. Shown, is a schematic of a seg- 
ment of the p**" row of Cu-O sites and its immediate neighbor- 
hood. Cu sites are located at the intersections of the lines; 
O sites are half-way between two Cu sites. Filled (empty) 
circles are occupied (empty) sites. Empty squares are sites 
that can either be occupied or empty. In a) we consider the 
change in the state produced when a hole hops between a Cu 
site labeled 1 in the row and the O site labeled 2. In the 
initial state, the Cu site must be occupied by a hole and the 
O site must be empty. In order that both the initial and fi- 
nal state survive projection with Pq, the O sites to the left 
and right must be empty and the O site just below must be 
occupied. Thus, this process necessarily decreases n{pax) by 
1, increases Nx{pax) by 1, and hence leaves Xp unchanged. 
In b) we illustrate the same considerations for applied to the 
state in which a hole hops from the O site labeled 3 to the Cu 
site labeled 4, and the notation is the same. The state of the 
O sites with squares is not uniquely determined. In this case, 
this process leaves unchanged n{pax) = 1 and Nx{pax) — 0, 
and so does not change Xp. Clearly, the same considerations 
imply that the inverse processes, and the processes involving 
sites immediately below the p"* row conserve Xp. All other 
processes in Hi trivially commute with Xp. 



ber of "quasiparticles" in each row or column: 

^P = Y1 Mi^p) + ^^il^P) + NA<1,P)} (A8) 



where {q,p) is the Cu site qax + pay, and Nx{R) is an 
operator defined by 

if ny{R)=n^y{R) = l, niR) = 



otherwise 



(A9) 
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Explicitly, since doubly occupied sites are anyway sup- 
pressed hy Ud and Up, 

N^iR) = [1 - niR)]ny{R)n^y{R). (AlO) 

The first two terms in the sum in Eq. l|A8p count the 
number of holes along the p*'' row. The third term, as 
we shall see, properly accounts for the finite transverse 
width of the actual quasi-particle excitations by counting 
the number of vacant Cu sites along the row which have 
an occupied O site both above and below. The column 
operators, Yp, and the associated operator Ny{R), are 
defined analogously. 

To show that each of these quantities is conserved, we 
first compute the commutator [Hi , Xp] and then (since 



[Po7-^p] = 0) sandwich the resulting expression between 
projection operators. Xp trivially commutes with all 
terms in Hi except those that hop a hole between a Cu 
site in row p and the O site immediately above (Fig. 
or below it, or between a neighboring O sites and the 
next Cu site immediately beyond it (Fig. ^.) In gen- 
eral, the application of these terms changes the value of 
Xp. However, as can be seen in the figure, and is ex- 
plained in the caption, if we enforce the condition that 
the initial and final state after the application of Hi is 
still a zero energy eigenstate of Hq, only processes which 
conserve Xp survive. This completes the proof. 
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